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Novel physical mechanism is proposed for explanation of temperature-independent transition re- 
actions in molecular systems. The mechanism becomes effective in the case of conformation transi- 
tions between qucisi-isoenergetic molecular states. It is shown that at room temperatures, stochastic 
broadening of molecular energy levels predominates the energy of low frequency vibrations accom- 
panying the transition. This leads to a cancellation of temperature dependence in the stochastically 
averaged rate constants. As an example, physical interpretation of temperature-independent onset 
of P2X3 receptor desensitization in neuronal membranes is provided. 
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Introduction.- Biomacromolecules (proteins and nu- 
cleic acids) operate by maintaining the set of conforma- 
tions and making transitions between them with proper 
reaction rates which rule, temperature con- 

trolled. Description of temperature-independent effects 
requires dissection of multi-type reactions arising within 
the complexity of molecular structures. In biology, one 
can generally indicate three types of reactions. The first 
type is generally associated with photoinduced oxidation- 
reduction processes Basic mechanism of these reactions 
is an electron tunneling between rigid redox centers like 
metal containing porphyrin rings and quinones The 
fact of tunneling is supported by the absence of temper- 
ature dependence for the transfer rates in wide tempera- 
ture regions including the room temperatures 0, Q . The 
second type of reactions is associated with enzyme kinet- 
ics [1] where the rate constants are calculated using the 
concept of activated complex. Motion of reactants along 
the reaction coordinates is accompanied by overcoming 
the activation barriers so that the rate constants experi- 
ence exponential temperature dependence. However, in 
the framework of Marcus's theory, a barrierless behav- 
ior becomes also possible if only reorganization energy 
and driving force of reaction are equal thus smoothing 
over temperature dependence. Such situation is proba- 
bly realized during the binding of NO molecule to the 
protoheme in the region of 200-290 K 5]. The third type 
of reactions is associated with conformation transforma- 
tions. Some of these reactions do not exhibit temperature 
dependence even in the range of room temperatures. For 
instance, it has been demonstrated that the closing rate 
for both the single-stranded DNA and RNA hairpin-loop 
fluctuations and the rate for the cyclic /? -hairpin pep- 
tide folding are weakly dependent on temperature in the 
interval from 10° to 60°C showing zero or even slightly 
negative enthalpies 0| . Another example of temper- 
ature independence in the gating kinetics of membrane 
proteins has been recently provided for the desenstization 
onset of P2X3 purinoreceptors [8] . 

While physical mechanisms explaining both 



temperature-independent tunnel electron transfer 
and barrierless ligand binding are more or less clear for 
two noted types of reactions, it is not the case for the 
third one. But, we have to note here a recent paper of 
Pouthier [9'] devoted to the study of exciton diffusion in a 
lattice of H-bonded peptide units. As it follows from the 
paper, a transition rate can become almost temperature 
independent in the strong anharmonic limit. Thus, the 
anharmonic nature of the phonons can be thought as one 
of possible mechanisms forming the specific temperature 
behavior of transition rates. 

In the present communication, we propose a 
novel physical mechanism explaining the existance of 
quasi-isoenergetic temperature-independent reactions in 
biosystems at room temperatures. It is assumed that for 
such type of reactions, the transitions between molec- 
ular states are accompanied by low-frequency molecu- 
lar vibrations, while molecular energies are alternated by 
thermodynamic fluctuations. 

Model and master equation.- Model is based on the ex- 
perimentally supported fact that protein-structural mo- 
tions create electric field fluctuations accompanying the 
movements of polar residues fl^. Such structure fluctu- 
ations cause the molecular energies to be stochastically 
time dependent and, thus, the transitions between molec- 
ular states occur on the background of random time- 
dependent energy shifts. Associating the molecule with 
an open quantum system S, the respective molecular 
Hamiltonian can be represented as H{t) — HQ{t) + V 
where iJo(i) — X^n Enit)\n){n\ is the main part of Hamil- 
tonian with En (t) = En + AEn (t) being the state energy 
varied owing to the stochastic addition A£'„(i). Transi- 
tions between molecular states \n) and \m) are achieved 
by nonadiabatic operator V. To analyze the experimen- 
tal results one have to know an evolution of observable 
occupancy Pn{t) for each molecular state involved into 
a given transition process. In our case, this occupancy 
appears as an average over realizations of random energy 
shifts. We denote such averaging via the symbol ((...}} so 
that Pn{t) = {{'Pn{t))). In turn, non-averaged occupancy 
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Vn{t) can be found from exact relation Pn{t) = {n\p{t)\n) 
where molecular density matrix p{t) — trBPs+B{t) is de- 
termined as a trace (over the thermal bath states) on 
density matrix ps+B{t) of an entire quantum system in- 
cluding, a thermal bath with respective bath Hamiltonian 
Hb [III ■ Evolution of the ps+ b {t) is governed by Liou- 
ville equation ps+Bit) = {l/h)[Ho{t) + V+HB, Ps+B{t)]. 
Using this equation one can derive an exact master 
equation just for the p{t). For our purposes, however, 
it is quite enough to have a master equation for only 
a diagonal part of molecular density matrix, Pd{t) = 
^„ (n|p(t)|n)|n)(n| = Dp{t). In what follow, we asso- 
ciate nuclear vibrations with a thermal bath. It allows 
one to use a well known factorization ps+B{t) = p{t)pB 
where ps — exp {—HB/kBT)/trB exp (—HB/kBT) is the 
bath equilibrium density matrix. Using now a method 
of projection operators 12] and generalizing this method 
for the case of time-dependent Hamiltonian, we achieve 
the following master equation 

Pd{t) = - dr 
Jq 

xtrB{D[V, U{t,t- t)[V, pd{t - T)pB]U+{t,t- r)]} (1) 

where the evolution operator reads U(t,t — t) = 
fexp[-{i/h)J^_^dT'[Ho{T') + HB + {l~D)V] (f is 
the Dayson's chronological operator). Eq. ([IJ is the basic 
coarse-grained stochastic equation which can be applied 
for a rigorous description of transitions between molecu- 
lar states on the time scale At ^ Tyu, where r„ib is the 
characteristic time of establishment of a thermal equilib- 
rium. 

Averaged transition rate constant- In what follows 
we concentrate our attention on the transition processes 
caused by a weak coupling of molecular states to the vi- 
brations. As it follows from a theory 1, 13], if energy of a 
single phonon is able to cover a difference between ener- 
gies of molecular states involved in the transfer process, 
the main contribution in transfer rate follows from sin- 
gle phonon-transitions while a role of multi-phonon pro- 
cesses becomes minor. Just such situation is assumed to 
be valid for low-energetic transitions in flexible molecules. 
The noted circumstances allow one to restrict a deriva- 
tion to the Born approximation over nonadiabatic oper- 

ator MV = En.m (1 - Km) Ex ^nm " ) |n) (m| . 

Here, (3\ ( (3~^) is the phonon annihilation (creation) op- 
erator for the Ath mode and = — is the cou- 
pling constant. Substituting the V in master equation ^ 
and taking the bath Hamiltonian in a convenient form, 
Hb = X^A + 1/2) {^\ is the frequency of Ath 

mode), one achieves the following set of non-Markovian 
equations 

m A 

x[e*""'""J-,„„(t,t-r)7'„(t-r) 
-e^"™"^ J-™„ {t,t - T)V^{t - t)] . (2) 



In this integro-diiferential equation, J7„„i = {En — 
Em)/fi is the transition frequency, J>„„(t,t — r) = 
exp{(i/?i) dT'[A£'„(T') - AE„,(t')\) is the stochas- 
tic functional, and R\{t) = N{uj\) exp {iuj\T) + [1 + 
N(u!x)] exp {—iui\T) is the regular single-phonon factor 
with N{ujx) = [exp{hw\/kBT) — being the phonon 
distribution function. 

To derive equations for the observable occupancies 
Pn{t), one has to average the stochastic equation ([2|). 
Here, we consider only the case of fast variations for each 
random energy difference A£'„(r) — AEm{T) when the 
characteristic times of stochastic and transition processes 
{Tstoch and Ttr, respectively) obey inequality Tstoch < Tjr- 
This inequality allows us to use a decoupling procedure 

{{Pnrn{t,t - T)Pn{t - t))) = {{Tnmit^t - t))) Pn{t - 

t). Note now that in the Born approximation, a non- 
Markovian character of the transition process becomes 
not important so that one can set P„ {t — T) « P„ [t) (see 
more details in ref. [l^). Estimation of averages like 
{{J^nm{t, t — t))) = Fnm{T) is givcu iu uumerous papers 
(see, e.g. refs. (l6l421| ). The results indicate that inde- 
pendently of a concrete form of quantities Fmn (t) , all of 
them decay with own characteristic times Tstoch- Taking 
this fact into account and using the condition Tstoch ^ Ttr 
one can shift the upper limit in the integral from t up to 
infinity. Thus, the averaging over the fast random real- 
izations transforms stochastic equation ([2]) into balance 
like kinetic equation, 

Pn (0 = - E [^^"™^" (*) - K,nnP,n {t)] (3) 

with the averaged transition rate constants 

= 1^ E \<m? dTe^^-^ Rx{T)Fn,n[T) . 

(4) 

Now we consider an important random process which 
allows for an exact calculation of the Fnm{T)- Let 
stochastically caused alternation of molecular energy dif- 
ferences be characterized by the only random quantity 
anm[t) so that A£;„(r) - AEmiT) = 27rto„m(T). For 
the sake of definiteness, let random quantity fluctu- 
ate between two equiprobable values +a and —a with 
mean frequency v (symmetric dichotomous process, DP) 
so that anm{i) = Oi{i) = if- Such random varia- 
tion of energy shifts yields Fnm[T) — F{t) — {{X{t))) 
where X(t) = exp [z JJ" dT'a{T')\ is the stochastic func- 
tional. Bearing in mind exact properties of the DP 
so as a^(r) — a^, ((q;(t))) = 0, and ((a(r)a(To))) = 
cr^ exp (— i^|t — ro|) and using the Shapiro-Loginov the- 
orem [17,] in the form of relation d{{a(T)X{T)))/dT ~ 
ivd{{X{T)))/dT + ia'^{{X{T))), one obtains F{t) = 
{kie~^'^'^ - k2e^^^'^)/{ki — ^2). Here, quantities ki^2 ~ 
v/2± ^ {v/^Y — cr^ determine the above noted charac- 
teristic stochastic times t^^'^^j^ = (Refci^2)^^- Substitut- 
ing the Fnm{T) = F{t) in Eq. Q and calculating the 
respective integral one achieves the following expression 
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for an averaged transition rate constant 

r 



+ {N{iJx) + l)L{-f,iy;n,„n-~ojx) ■ (5) 



Here, the influence of stochastic field on the n ^ m 
transition is concentrated in Lorentzian like function 
L{'y,i/;n) — TT^^j/[{'y — ft^/h')'^ + fl^] via parameters 
u and 7 = (T^/i/. 

Temperature independence of rate constant— Rate con- 
stant ([5]) includes a stochastic field in a non-perturbation 
manner. This circumstance allows us to analyze various 
regimes of transition processes in molecular systems de- 
pendently on value of field parameters v and 7. Here, 
bearing in mind an application of theory to explanation of 
conformational quasi-isoenergetic transitions (fi„m ~ 0) 
in fiexible biomolecules at room temperatures, we re- 
strict ourself to the analysis of transitions which are ac- 
companied by low frequency vibration modes so that 
\flnm i 1^- This rcduccs function L{'y,i';fl) to 

a standard Lorentzian L(7, ft) « 7r~^7/(7^ -I- ri^) where 
parameter /17 can be treated as the broadening of molecu- 
lar energy levels. At 7 — )■ 0, i.e. in the case of extremely 
high frequency stochastic field {v — ?> 00), Lorentzian is 
converted into the Dirac's delta-function S{fl) and, thus, 
rate constant ([5]) appears in a conventional form valid for 
Born approximation over a nonadiabatic perturbation. 
In the case of moderately high frequency stochastic field 
under consideration, an another limit is realized where 
7 ^ ri. It yields ^(7, ^) = {ttj)~^. Note now that room 
temperatures correspond to energies of the order 0.025 
eV, identical to 6 • lO^^s^^ or 200 cm^^. Therefore, if 
transitions are accompanied by the vibrations of the or- 
der 60 cm~^ and lower, one can set N{u}x) ~ kBT/hux- 
This reduces averaged rate constant ([5]) to the form 
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(6) 



Physical origin of stochastic parameter 7 is dictated 
by random shifts of molecular energy levels. We pro- 
pose a phenomenological model where mean positive and 
mean negative shifts (+cr/2 and — cr/2, respectively) re- 
sult from thermodynamic fluctuations of an energy so 
that (cr/2)^ — SE"^ /{2TThY . In line with a general theory 
of thermodynamic fluctuations in canonic ensembles, the 
aver age of square of energy fluctuations is calculated as 
= ksT'^idE/dT) [22] with 'E being the mean en- 
ergy of a system. At the same time, the average fre- 
quency of fluctuations is associated with a linear fre- 
quency V — E/2TTh \4^. Thus, a stochastic parameter can 
be estimated with relation 7 = {2kBT^ /TTh){dlnE/dT). 
In classical limit under consideration, a mean energy per 
a separate degree of freedom, is exactly equal to thermal 
energy /c^T. Therefore, independently of precise molec- 
ular system, a total mean energy E is proportional to 
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FIG. 1: Temperature independence of desensitization onset of 
P2X3 receptors measured at physiological temperatures 25° C- 
40° C (selected from @]) and its fit by Eq. ([8]). Insets indicate 
limiting cases in the model of two types of conducting chan- 
nels (details in the text). 



a bath temperature, 
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It yields 7 — 2kBT/nh. Conse- 
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(7) 



In Eq. dZl), i?„™ = {2TThy^ J^^ {duj/uj) Jr,m{uj) is the 
temperature-independent factor which is expressed via 



spectral function Jnm(s^) — '^'^^x 



'5(wA-w). The 



latter is widely used in a spin-boson model for descrip- 
tion of transport processes accompanied by vibration mo- 
tions . In a given case, function J„m(a;) includes 
an information on both a vibrational structure of fiex- 
ible molecule and a character of nonadiabatic coupling 
to molecular vibrations. 

Experimental evidence.~Eq. ([7]) contains a funda- 
mentally important result indicating that the rates of 
quasi-isoenergetic transition processes in molecules can 
be temperature- independent even at room temperatures. 
As an example, we consider the desensitization onset 
of P2X3 receptors in neuronal membranes [8]. For 
our knowledge, this is the first quantitative observation 
of temperature-independent gating in biological mem- 
branes. Fig. [1] manifests the decrease of ion current I{t) 
through the receptor channels. Desensitization probabil- 
ity of the channels, Pd{t) — 1 — /(t)//(0) is well described 
by a two-exponential kinetics Q, 



Pd{t) = l-Aie-'/''~A2e-'^^^ 



(8) 



where ti ~ 14.7ms and T2 ~ 231ms are the temperature- 
independent characteristic times of desensitization while 
Ai ~ 0.968 and A2 ~ 0.032 are the pre-exponential 
weights. Since a particular molecular content of the P2X3 
receptor gates is yet unknown, we interpret the noted re- 
sult in the framework of simplest model supposing the 
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formation of ion current as the sum of two partial cur- 
rents. Let ij be the current through a separate open 
channel of the j{— l,2)th type. If Nj is the number of 
respective ion channels, then a partial current reads as 
Ij{t) = Nj ijP^^\t) where Py\t) is the probability for 
the channel to be in the open state "o". Desensitiza- 
tion of each channel occurs independently of one another 
manifesting the process of conformational transition from 
the open (conductive) state into the closed (nonconduc- 
tive) state so that Po^\t) = exp{—t/Tj). Introducing 
the quantities /(O) = J2j=i 2 ^3 ^'^'^ ^1 — ^/(l + 
and A2 = 1/(1 + ^) where ^ = A^i ii /A^2 ^2 , one ex- 
presses an ion current in the form I{t) — 1(0) Po{t) where 
Po{t) = X]j=i 2 Aj pjy-'\t) is the apparent (statistically 
averaged) probability of an ion channel to be in the open 
state. Bearing in mind that Pd{t) = I — Po{t), one arrives 
at the Eq. ©. 

Physical explanation of temperature independency of 
desensitization process can be given on the base of above 
proposed model of quasi-isoenergetic transitions in flex- 
ible molecules. Denoting via n = joa and m = jda' 
the conformational isoenergetic substates participating 
in the open-close transition, we introduce the integral 

occupancies of states o and d as Po''\t) = Ylia=i Pjoa{t) 
and P^'\t) = Y.T=i Pidc'it) = 1 - P^o'\t) where /x^o 
and /ijd are the number of respective degenerated sub- 
states for the jth type of channel. Eq. ^ reads 

now as P^o'\t) = - [Ki'ldP^o'\t) ~ K^dloPd\t)] ^here 
^l^id — Kj/^jo and K^\^ = Kj/^jd are the integral 
rate constants with ivTj = {2Tr/h) J2a=i Sa'li Rjoajda'- 
Since the disordering increases the degeneracy of molec- 
ular state, then pLjd ^ /ijo- This reduces the ki- 
netics of isoenergetic transitions to the nonrecurrent 
single-exponential kinetics with the characteristic time 



Tj ~ {Kol,d)~^^ fo'" the jth type of channel. Thus, 
physically, the receptor desensitization appears as a 
temperature-independent open— i>close transition process 
between quasi-isoenergetic molecular conformations of 
ion channel. 

Conclusion.-ln this communication, a novel physi- 
cal mechanism is proposed to explain the formation of 
temperature-independent transition processes in molec- 
ular systems. In contrast with well established mecha- 
nism describing quantum phononless site-to-site particle 
tunneling at low temperatures or tunneling at room tem- 
peratures that is accompanied by an emission of high 
frequency phonons (when wa ^ fcBr/ft)[l[, the pro- 
posed mechanism works in a classic region of tempera- 
tures with participation of low frequency phonons (when 
i^x ^ ksT/h). The lack of temperature dependence in 
rate constant ^ occurs owing to the thermodynamical 
stochastic variation of energy levels participating in the 
transitions. Frequency of this variation, v is assumed 
to be fast in comparison to the transition frequencies 
flnm ^ This reduces an averaged rate constant ([5]) 
to a more simple form ^ which, in turn, is simplified 
to the temperature-independent form ([7]). Of course, a 
given mechanism is not unique and can work indepen- 
dently of or in parallel with other possible mechanisms, 
especially with anharmonic mechanism noted in ref. @. 
But, for understanding the physics of quasi-isoenergetic 
transitions in biological macromolecules, the reference of 
fast molecular motions to the motions responsible for the 
creation of high frequency stochastic fields, is assumed to 
be quite fruitfuU from a semi-phenomenological physical 
point of view. 
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